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Abstract 

Different perturbation theory treatments of the Ginzburg-Landau phase transition 
model are discussed. This includes a criticism of the perturbative renormalization 
£h ' group (RG) approach and a proposal of a novel method providing critical exponents 

consistent with the known exact solutions in two dimensions. The new values of 
critical exponents are discussed and compared to the results of numerical simulations 
and experiments. 
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O ' 1 Introduction 

O 

Phase transitions and critical phenomena is one of the most widely investigated top- 
ics in modern physics. Nevertheless, a limited number of exact and rigorous results is 
available Our purpose is to give a critical analysis of the conventional approach in 
calculation of critical exponents based on the perturbative renormalization group (RG) 
theory ||, [3|, ||] and to propose a new method which provides results consistent with the 
known exact solutions. The basic hypothesis of the conventional (RG) theory is the ex- 
istence of a certain fixed point for the RG transformation. However, the existence of 
such a stable fixed point for the Ginzburg-Landau model (which lies in the basis of the 
^ ■ field theory) has not been proven mathematically in the case of the spatial dimensionality 

d<4. 

The usual RG theory treatment of the Ginzburg-Landau model is based on the di- 
agrammatic perturbation theory (Feynman diagrams). We have demonstrated that this 
treatment is contradictory and therefore cannot give correct values of critical exponents. 
Namely, based on a method which is mathematically correct and well justified in view of 
the conventional RG theory, we prove the nonexistence of the non-Gaussian fixed point 
predicted by this theory (Sect. |||). In Sect. || we prove that a correctly treated diagram 
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expansion provides results which essentially differ from those of the perturbative (diagram- 
matic) RG theory. Finally, we have proposed a novel analytical method of determination 
of critical exponents in the Ginzburg-Landau model (Sect. ||), and have compared the pre- 
dicted exact values of critical exponents to the results of numerical and real experiments 
(Sect. D- 

2 Critical analysis of the perturbative RG method 

Here we consider the Ginzburg-Landau phase transition model within the usual renor- 
malization group approach to show that this approach is contradictory (for more details 
see also [ 10 1 ) . The Hamiltonian of this model in the Fourier representation reads 



T =H( r ° + ck2 ) I \ 2 + uV 1 ^kx^ka^a^-ki-ka-ka > i 1 ) 

k ki,k 2 ,k 3 

where ip^ = V~ l l 2 J p(x) exp(— zkx) cbc are Fourier components of the scalar order pa- 
rameter field </?(x), T is the temperature, and V is the volume of the system. In the RG 
field theory ||, ||| Hamiltonian (|l|) is renormalized by integration of exp(—H/T) over tp^ 
with A/s < k < A, followed by a certain rescaling procedure providing a Hamiltonian 
corresponding to the initial values of V and A, where A is the upper cutoff of the ip 4 
interaction. Due to this procedure, additional terms appear in the Hamiltonian (|l|), so 
that in general the renormalized Hamiltonian contains a continuum of parameters. The 
basic hypothesis of the RG theory in d < 4 dimensions is the existence of a non-Gaussian 
fixed point \i = /i* for the RG transformation R s defined in the space of Hamiltonian 
parameters, i.e., 

RsV* = M* • (2) 

The fixed-point values of the Hamiltonian parameters are marked by an asterisk (rg, c* , 
and u*, in particular). Note that /i* is unambiguously defined by fixing the values of c* 
and A. According to the RG theory, the main terms in the renormalized Hamiltonian in 
d = 4 — e dimensions are those contained in (jl|) with 7q and u* of the order e, whereas the 
additional terms are small corrections of order e 2 . 

Consider the Fourier transform G(k, ^u) of the two-point correlation (Green's) function, 
corresponding to a point /i. Under the RG transformation R s this function transforms as 
follows § 

G(k, fj.) = s 2 ^ G(sk, R s n) . (3) 

Let G(k, /i) = G(k, fi) (at k ^ and V — > oo) be defined within k < A. Since Eq. (||) 
holds for any s > 1, we can set s = A/k, which at /x = /x* yields 

G(k, /j*) = a k- 2+r > for k < A , (4) 

where a = A 2 ~ V G(A, /i*) is the amplitude and rj is the universal critical exponent. Ac- 
cording to the universality hypothesis, the infrared behavior of the Green's function is 
described by the same universal value of rj at any /i on the critical surface (with the only 
requirement that all parameters of Hamiltonian (|l]) are present), i.e., 

G(k, fi) = b(n) k~ 2+v at k -» , (5) 
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where 

%) = hm^G(k, / u). (6) 

re— >0 

According to Eq. (||), which holds for any s = s(k) > 1 and for s = A/k in particular, 
Eq. (|6j) reduces to 

b(ji) = lim k^sikf-* G(sk, R s fi) = a , (7) 
if the fixed point u* = lim R s fi exists. Let us define the function X(k, u) as A(k, fi) = 

s — >oo 

k~ 2 G~ 1 (k, n). According to Eqs. (||), @, and (0), we have (for A; < A) 

X(k,//) = -AT^ (8) 
a 

and 

X(k,/z) = -AT" + 5X(k,/x) , (9) 
a 

where /x belongs to the critical surface, and <5X(k, fi) denotes the correction-to-scaling 
term. From (||) and (|9|) we obtain the equation 

5X(k, n* + dfi) = X(k, /x* + 5fi) - X(k, /x*) , (10) 

where <5/x = /x — /x*. Since Eq. ( |l0| ) is true for any small deviation 5/tx satisfying the relation 

/x* = lim R s (fi* +6fj.) , (11) 

S — >OG 

we choose <5/i such that /j,* ^ n* + b[i corresponds to the variation of the Hamiltonian 
parameters 7q =^ r$ + <5ro, c* =^ c* + <5c, and u* 4«* + 6xA, where A is a small constant. 
The values of 5r$ and 5c are choosen to fit the critical surface and to meet the condition 
(|PI|) at fixed c* = 1 and A = 1. In particular, quantity 5c is found 5c = B e 2 + o(e 3 ) with 
some (small) coefficient B = i?(A), to compensate the shift in c of the order e 2 due to the 
renormalization (cf. ||]). The formal e-expansion of 5X(k, /x), defined by Eq. (|l0|), can be 
obtained in the usual way from the perturbation theory. This yields 

5X(k,fx) = e 2 [Ci(A) +C 2 (A) lnfc] + o(e 3 ) at fe -> , (12) 

where Ci(A) and C 2 (A) (C 2 7^ 0) are coefficients independent on e. 



It is commonly accepted in the RG field theory to make an expansion like (12), obtained 
from the diagrammatic perturbation theory, to fit an asymptotic expansion in k powers, 
thus determining the critical exponents. In general, such a method is not rigorous since, 
obviously, there exist such functions which do not contribute to the asymptotic expansion 
in k powers at k — * 0, but give a contribution to the formal e-expansion at any fixed 
k. Besides, the expansion coefficients do not vanish at k — > 0. Trivial examples of such 
functions are e m exp(— ek~ e ) and e m [1 — tanh(efc -e )] where m is integer. Nevertheless, 
according to the general ideas of the RG theory (not based on Eq. (10)), in the vicinity 
of the fixed point the asymptotic expansion 



X(k, n) = -k-v + hk e+ °^ + b 2 k 2+ °^ + ... (13) 
a 

is valid not only at k — > 0, but within k < A. The latter means that terms of the kind 
e m exp(—ek~ e ) are absent or negligible. Thus, if the fixed point does exist, then we can 
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obtain correct e-expansion of 5X(k,fi) at small k by expanding the term bik e+ °^ (with 
b\ = b\ (e, A)) in Eq. (|l3| ) in e powers, and the result must agree with (|l2| ) at small A, at 
least. The latter, however, is impossible since Eq. (12) never agree with 



<Hf(k,/i) =&!(€, A) 



1 + elnk + o(e" 



(14) 



obtained from ( [131 ) at A; — » 0. Thus, in its very basics the perturbative RG method in 4 — e 
dimensions is contradictory. From this we can conclude that the initial assumption about 
existence of a certain fixed point, predicted by the RG field theory in 4 — e dimensions, is 
not valid. 



3 A model with quenched randomness 

Here we consider the Ginzburg-Landau phase transition model with 0(n) symmetry (i.e., 
the n-vector model) which includes a quenched randomness, i.e., a random temperature 
disorder (for more details see also |11[| ). One of the basic ideas of the perturbative RG 
theory is that n may be considered as a continuous parameter and the limit n — > 
makes sense describing the self-avoiding random walk or statistics of polymers [Q, ||]. We 
have proven rigorously that within the diagrammatic perturbation theory the quenched 
randomness does not change the critical exponents at n — > 0, which is in contrast to the 
prediction of the conventional RG theory formulated by means of the Feynman diagrams. 
The Hamiltonian of the actually considered model is 

H/T = |[(r + ^/(x))^ 2 (x) + C (V 95 (x)) 2 ]dx (15) 

+ uV' 1 <M k i)</>i( k 2) w kl +k 2 <Pj(k 3 )tpj(-ki - k 2 - k 3 ) 

i.j,ki,k 2 ,k 3 

which includes a random temperature (or random mass) disorder represented by the term 
y/u /(x) y 2 (x). For convenience, we call this model the random model. In Eq. ( JT5| ) y(x) 
is an n-component vector with components </?i(x) = V~ 1 / 2 J2k<A Vi( k ) e * kx > depending 
on the coordinate x, and /(x) = V~ 1 / 2 X)k /ke* kx is a random variable with the Fourier 
components /k = V^ 1 / 2 J /(x)e~* kx dx. 

The system is characterized by the two-point correlation function Gj(k) defined by 
the equation 

(^(k)^.(-k)) = tfy Gi(k) = 5 itj G(k) . (16) 

It is supposed that the averaging is performed over the </j( x ) configurations and then over 
the /(x) configurations with a fixed (quenched) Gaussian distribution P({/k}) for the set 
of Fourier components {/k}; i- e., our random model describes a quenched randomness. 
We have proven the following theorem. 

Theorem. In the limit n — ► 0, the perturbation expansion of the correlation function 
G(k) in u power series for the random model with the Hamiltonian ( |15"D is identical to the 
perturbation expansion for the corresponding model with the Hamiltonian 

H/T = J [ro^ 2 (x) + c(V^(x)) 2 ] dx (17) 

+ uV' 1 Y ^( k i)^i( k 2) ^ki+k 2 Pj(k3)¥>j(-ki - k 2 - k 3) 

i,j,ki,k 2 ,k 3 
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where u k = u k - \ ^| / k 

For convenience, we call the model without the term yju /(x) 9? 2 (x) the pure model, 
since this term simulates the effect of random impurities || . 

Proof of the theorem. According to the rules of the diagram technique, the formal 
expansion for G(k) involves all connected diagrams with two fixed outer solid lines. In the 
case of the pure model, diagrams are constructed of the vertices ^>«m**<^ , with factor 
— ny _1 nk related to any zigzag line with wave vector k. The solid lines are related to the 
correlation function in the Gaussian approximation Go(k) = 1/ (2ro + 2ck 2 ). Summation 
over the components ifi (k) of the vector V'(k) yields factor n corresponding to each closed 
loop of solid lines in the diagrams. According to this, the formal perturbation expansion 
is defined at arbitrary n. In the limit n — > 0, all diagrams of G(k) vanish except those 
which do not contain the closed loops. In such a way, for the pure model we obtain the 
expansion 

G(k) = k ± + krA-k + ... . (18) 



In the case of the random model, the diagrams are constructed of the vertices 
and ••••< . The factors uV 1 / k | 2 ^> correspond to the coupled dotted lines and the 
factors -mF" 1 ^ correspond to the dashed lines. Thus, we have 



G(k) = k :k + 



+ ... . (19) 



In the random model, first the correlation function G(k) is calculated at a fixed {/ k } 
(which corresponds to connected diagrams where solid lines are coupled, but the dotted 
lines with factors —y/uV^^f-^ are not coupled), performing the averaging with the weight 
-P({/k}) over the configurations of the random variable (i.e., the coupling of the dotted 
lines) afterwards. According to this procedure, the diagrams of the random model in 

general (not only at n — > 0) do not contain parts like , ^ , , etc., 

which would appear only if unconnected (i.e., consisting of separate parts) diagrams would 
be considered before the coupling of dotted lines. 

It is evident from Eqs. ( |l8| ) and ( |l9| ) that all diagrams of the random model are 
obtained from those of the pure model if any of the zigzag lines is replaced either by 
a dashed or by a dotted line, performing summation over all such possibilities. Such a 
method is valid in the limit n — > 0, but not in general. The problem is that, except the 

case n — > 0, the diagrams of the pure model contain parts like @$r~* , ^CZ. > ^6cE > 
etc. If all the depicted here zigzag lines are replaced by the dotted lines, then we obtain 
diagrams which are not allowed in the random model, as explained before. At n — > 0, 
the only problem is to determine the combinatorial factors for the diagrams obtained by 

the above replacements. For a diagram constructed of M\ vertices ^> <C and M2 

vertices the combinatorial factor is the number of possible different couplings of 

lines, corresponding to the given topological picture, divided by M\\M<2}.. 

Our further consideration is valid also for the diagrams of free energy (at n — > 
represented by the main terms containing single loop of solid lines) and of 2m-point 
correlation function. We define that all diagrams which can be obtained from the z-th 
diagram (i.e., the diagram of the i-th topology) of the pure model, belong to the i-th. 
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group. Obviously, all diagrams of the i-th group represent a contribution of order u l , 
where I is the total number of vertices >^^< in the i-th diagram. The sum of the 
diagrams of the i-th group can be found by the following algorithm. 

1. Depict the i-th diagram of pure model in an a priori defined way. 

2. Choose any one replacement of the vertices >*ww<^ by > <C and ]>• *C' 

and perform the summation over all such possibilities. For any specific choice we 
consider only one of the equivalent M\\M%\ distributions of the numbered M\ ver- 
tices >" <C and M2 vertices ••••< over the fixed numbered positions instead 

of the summation over all these distributions with the weight \j{M.]\M.<£). Thus, 
at this step the combinatorial factor for any specific diagram is determined as the 
number of possible distributions of lines (numbered before coupling) for one fixed 
location of vertices consistent with the picture defined in step 1. 

3. The result of summation in step 2 is divided by the number of independent symmetry 
transformations (including the identical transformation) for the considered i-th dia- 
gram constructed of vertices ^>w~\<^ , since the same (original and transformed) 
diagrams were counted as different. 

Note that the location of any vertex >» <C is defined by fixing the position of dashed 

line, the orientation of which is not fixed. According to this, the summation over all 
possible distributions of lines (numbered before coupling) for one fixed location of vertices 
yields factor 8 Ml 4 Ma / 2 . The i-th dia gram of the pure model also can be calculated by 
such an algorithm. In this case we have 8 ; line distributions, where I = M\ + M2/2 is 
the total number of vertices ]>*w«^<; in the i-th diagram. Obviously, the summation 
of diagrams of the i-th group can be performed with factors 8' instead of 8 Ml 4 M2 / 2 , 
but in this case twice smaller factors must be related to the coupled dotted lines. The 
summation over all possibilities where zigzag lines are replaced by dashed lines with factors 
— uy -1 Uk and by dotted lines with factors ^uV^ 1 (\ /k | 2 ^}, obviously, yields a factor 

uV~ l f— Uk + I (\ /k = — uV~ l U]s_ corresponding to each zigzag line with wave vector 
k. Thus, the sum over the diagrams of the i-th group is identical to the i-th diagram 
of the pure model defined by Eq. (|l7|). By this the theorem has proved not only for the 
two-point correlation function, but also for 2m-point correlation function and free energy. 

If, in general, the factor y/u in Eq. ([D|) is replaced by Vv/, where v! is an independent 
expansion parameter, then our analysis leads to the above relation between diagrams for 
utilt = uiik — y (\ /k | 2 ^- According to this, at n — ► the pure and random models 
cannot be distinguished within the diagrammatic perturbation theory. If, in principle, 
critical exponents can be determined from the diagram expansions at n — > 0, as it is 
suggested in the usual RG theory, then the same critical exponents should be provided 
for both models at n — > 0. In such a way, we conclude that the RG method is not 
correct because the above condition is violated. As compared to our simple treatment 
of the random model, the RG treatment includes additional Feynman diagrams because 
the Hamiltonian becomes more complicated after the renormalization. However, this does 
not enable to find the difference between both models: the original information, when one 



starts the perturbative renormalization of Hamiltonian (15), is contained in the Feynman 



diagrams we considered, but the renormalization by itself does not create new information 
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about the model. Really, by renormalization we merely "forget" some information about 
the short-wave fluctuations to make that for the long-wave fluctuations easier accessible. 
Thus, our conclusion remains true. 



4 New method based on perturbation theory 

As we have already discussed in Sect. it is not a rigorous method to make a formal 
expansion like ( |l2| ) and to try calculate the critical exponents therefrom. We propose 
another treatment of the diagrammatic perturbation theory. The basic idea is to obtain 
suitable equations by appropriate grouping of the diagrams. Suitable are such equations 
which allow to find the asymptotic expansions at the critical point directly in k power 
series, but not in terms of the formal parameter In k (as in Eq. |i~2] ) which diverges at k — ► 0. 
In such a way, for the Ginzburg-Landau model defined by Eq. (|i~7|), where jj% = u^, (for 
simplicity here we consider the case of the scalar order parameter, i. e. n = 1) we have 
obtained the Dyson equation 

1 2 dD(G) , s 

- r + ck 2 - ) ' (20) 



2G(k) KJ dG{k) 

where D{G) denotes a quantity, the diagram expansion of which involves all the so- 
called skeleton diagrams (constructed of the fourth-order vertices ^> <C with factors 

— y -1 ^ related to the dashed lines) without outer lines. Skeleton diagram is defined as 

a connected diagram containing no parts like ^ . In distinction to the usual 

(simple) perturbation theory, the true correlation function G(k) is related to the coupled 
solid lines instead of the Gaussian correlation function Go(k). In such a way, quantity 
D(G) have to be considered as a function of discrete variables G(k) corresponding to the 
set of discrete wave vectors k. Eq. fl20| ) has been previously obtained in Ref. |J|, where also 
the simplest (ring) skeleton diagrams have been considered. We have found a possibility to 
include all the skeleton diagrams into consideration, which allows to find the set of possible 
values for exact critical exponents. At u — > our equations, represented by converging 
sums and integrals, define the true correlation function G(k) with an error smaller than 
u l at any positive I. This is quite enough to find the exact critical exponents based on 
general scaling properties of the solution in vicinity of the critical point. The structure 
of our equations allows to prove these properties taking into account all diagrams. The 
results also can be easily generalized to the case with O(n) symmetry. It is not possible 
in this relatively short paper to give the mathematical derivations of our equations and 
results. Hopefully, they will be available in nearest future. Here we present the final result 
according to which possible exact values of critical exponents 7 (susceptibility exponent) 
and v (correlation length exponent) for the n-component vector model (n = 1, 2, 3, etc.) 
in d < 4 dimensions (i. e., at d = 2, 3) are given by 

= d+2j + 4m = 2(l+m)+j 

7 d(l + m + j) - 2j ' U d(l + m + j) - 2j ' 1 ' 

where m may have a natural value starting with 1 and j is integer equal or larger than 
— m. In general, different values of j and m can correspond to different (natural) n, i. e., 
j = j(n) and m = m(n). It is easy to verify that at j = and m = 3 Eq. (|2l| ) reproduces 
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the known Q exact results in two dimensions. The known exact exponents for the spherical 
model (n = oo) are obtained at j(n)/m(n) — > oo. Although the derivations are true for 
d < 4, Eq. provides correct result v = 1/2 and 7 = 1 also at d = 4. It is reasonable to 
consider (i as a continuous parameter. This leads to the conclusion that m = 3 and j = 
are the correct values for the case n = 1 not only at <i = 2, but also at <i = 3. In the latter 
case we have 7 = 5/4 and v = 2/3. The nearest values of 7 and v provided by Eq. (21), 
e. g., at j = 1 and m = 3 or at j = 1 and m = 4 are then the most probable candidates 
for the case n = 2. 



5 Comparison of results and discussion 

It is commonly believed that all more or less correct Monte Carlo (MC) simulations confirm 
the values of critical exponents obtained from the perturbation expansions based on the 
renormalization group. This is not true. We have found that some kind of MC simulations 
at the critical point, namely, the MC simulations of fractal configurations of Ising model || 
and the MC simulations of the energy density [[?J for the XY model in reality do not confirm 
the results of the RG theory, but provide the values of critical exponents which are very 
close to those we predicted. 

The MC simulations of Ref. [|| allows to determine the fractal dimensionality D (the 
largest cluster in the relevant configuration has the volume L D where L denotes the linear 
size of the system) which is related to the critical exponents by 7 = v(2D — d) or, which 
is the same, 77 = 2 — ~f /v = d + 2 — 2D. In our opinion, this method is better than other 
more convenient simulation methods, since it provides the value of rj as a result of direct 
simulation, i. e., there are no fitting parameters. Besides, the result is relatively insensitive 
to the precise value of the critical coupling (temperature). In Fig. [I] we have shown the 
average values of D (the averaging is is made over the MC steps from 1 to 10 (except 
the initial point), from 11 to 20, and so on) calculated from the MC data of Ref. || by 
measuring deviation from the line D = 2.48 in Fig. 8 (of Ref. ||). If properly treated, 
these simulation data confirm the value of r\ about 1/8 (or D = 2.4375) consistent with 
our prediction 7 = 5/4 and v = 2/3, as it is evident from Fig. |l|. The value D = 2.46 ±0.01 
reported in Ref. Q seems to be determined from the upper MC points (Fig. 8 in Ref. ||) 
only which are closer to the known theoretical prediction D = 2.48. 

As regards the MC simulations of the energy density E of XY model [Q] at the critical 
point, the true picture can be reconstructed from the simulated values listed in Tab. I 
of Ref. |0]. Since all the values of E are of comparable accuracy, it is purposeful to use 
the least-square method to find the optimum value of \/v by fitting the MC data to the 
prediction of the finite-size scaling theory 

E{L) = E + E x L-v~ d , (22) 

where E(L) is the energy density at the critical temperature T\ depending on the linear size 
of the system L. The standard deviation of the simulated data points from the analytical 
curve ( p2|) can be easily calculated for any given value of \jv with the parameters -Eo 
and E\ corresponding to the least-square fit. The result is shown in Fig. |2[ The thick 
solid curve is calculated including all 11 data points (L=10, 15, 20, 25, 30, 35, 40, 45, 
50, 60, 80), whereas the dashed line - including 9 data points (except L=10, 15) used 
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Figure 1: Fractal dimensionality D of the three dimensional Ising model at the critical 
point simulated by Monte Carlo method (MCS means Monte Carlo steps). The upper and 
lower dashed lines indicate the theoretical values expected from the known and from our 
critical exponents, respectively. 
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Figure 2: The standard deviation a vs the value of \jv used in the least-square fit of the 
finite-size scaling curve to the simulated results including 11 data points (solid curve) and 
9 data points (dashed curve). Minimum of the solid curve, shown by a vertical dashed 
line, corresponds to the best fit 1/v = 1.4457 which is close to our theoretical value 13/9 
indicated by a vertical dotted line. Other vertical dashed line indicates the value 1.487 
proposed by authors of Ref. |7| . 
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Figure 3: Our fits to the original MC data of Ref. @ for the energy density (left) and for 
the specific heat (right) depending on the linear size of the system L. 



for the fitting in Ref. @. Minimum of the solid curve, shown by a vertical dashed line, 
corresponds to the best fit \jv = 1.4457 which comes very close to our theoretical value 
13/9 ( provided by ( |2l| ) at j = 1 and m = 3) indicated by a vertical dotted line. We have 
estimated the statistical error of this MC result about ±0.007 by comparing the best fits 
for several random data sets. Different data sets have been generated from the original 
one by omitting some data points with 10 < L < 80. We have found it unreasonable to 
omit the data points with two smaller sizes, as it has been proposed in Ref. |7|, since the 
result in this case becomes very poorly defined, i. e., the dashed curve in Fig. | has a very 
broad minimum. Besides, there is no reason to omit the smallest sizes, since the analytical 
curve (|22[ ) excellently fit all the data points and the standard deviation for 11 data points 
is even smaller than that for 9 data points (see Fig. ||). The possible systematical error 
due to the inaccuracy in the critical temperature T\ = 2.2017 ± 0.0005 (the error bars are 
taken from the source of this estimation |||) used in the simulations |?J has been evaluated 
±0.017 by comparing the simulation results at T\ values 2.2012, 2.2017, and 2.2022. In 
this case the values of the energy density at a slightly shifted temperature have been 
calculated from the specific heat data given in Tab. I of Ref. j/J. In such a way, our final 
estimate from the original MC data of Ref. 0] is 1/v = 1.446 ± 0.025 in a good agreement 
with our theoretical value 13/9 = 1.444... and in a clear disagreement with the usual (RG) 
prediction about 1.493. One can only wonder where the value 1.487 proposed in Ref. 
comes from. It does not correspond neither to the best fit for 11 data points nor to that 
for 9 data points, as it is evident from Fig. ||[ The values of \jv and a/v cannot be 
determined independently from the discussed here energy density data. One of them have 
to be calculated from the scaling relation ajv + d = Ijv. If authors of Ref. were able 
to determine \ jv with ±0.081 accuracy, then they should be able to find a/v with ±0.162 
accuracy. In this aspect, the estimate a/v = —0.0258 ± 0.0075 given by the authors looks 
more than strange. 

In Fig. I we have shown our fits to the MC data for the energy density E(L) = 
2.0108 - 2.0286 L -14 / 9 and for the specific heat c(L) = 7.360 - 6.990 L" 1 / 9 . They do not 
look worse than those in Ref. j/J, but our fit for c(L) seems to be better. 

One believes that the value of critical exponent v about 0.67, predicted by the RG 
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Figure 4: Percent deviation of the experimental p s /p data || from the expected theoretical 
relation (Eq. (|23| ) at 5 = 0). The stright line shows the slope of this plot at the value of 
t, equal to 5 • 10 -7 , indicated by a vertical dashed line. 



theory at n = 2, is well confirmed by very accurate measurements of the superfluid fraction 
Ps/p = V in A He. This is not true, since in reality these experiments || provide a good 
evidence that the effective critical exponent v e ff(t) = d(lny)/d(lnt) remarkably increases 
when the reduced temperature t = (T\ — T)/T\ (where T\ is the critical temperature) is 
decreased below 10~ 5 . According to Ref. ||, p s /p is given by 

Ps/p = yit) = ko(l + ht)(l + D p t A )t< x (1 + 6(t)) , (23) 

where ko, ki, D p , and C, are the fitting parameters, A = 0.5 is supposed to be the 
correction-to scaling exponent, and 5{t) is the measured relative deviation from the ex- 
pected theoretical expression obtained by setting 5(t) = 0. The percent deviation discussed 
in Ref. ||] is 100 times S(t). From Eq. (p3|) we obtain 

. . ht AD p t A 1 dSCt) . t . 



For the values of t as small as t < 10 5 and for 6(t) <C 1 Eq. (p4]) with the fitting 
parameters £ = 0.6705, k$ = 2.38, k% = —1.74, and D p = 0.396 used in Ref. |]] reduces to 

u eff (t)~C + d6(t)/d(hit). (25) 

The second term in this equation is proportional to the slope of the percent deviation 
plot 100 5(t) vs Int or Igt (the decimal logarithm) in Figs. 2 and 3 of Ref. ||]. We have 
read the experimental data from Fig. 2 in Ref. [|| within the region t < 10 -4 and have 
depicted them in Fig. |I[ Almost all the data points with a reasonable accuracy fit the 
smooth curve 5(t) vs lg £ (dashed line) having a maximum at about lgi = —5.5. It means 
that d 2 5(t)/d(lnt) 2 is negative within some region around the maximum, i. e., according 
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to (pq) the effective critical exponent u e ff(t) increases if t is decreased. We have roughly 
estimated and have shown by stright line the slope of this curve at t = t* = 5-1CT 7 (t* value 
is indicated in Fig. || by vertical dashed line). From this we obtain d5(t)/d(lnt) ~ 0.025. 
This result depends on the shift in the experimentally determined T\ value. To obtain a 
more reliable estimate, we have performed the same manipulations with the data depicted 
in Fig. 3 of Ref. M corresponding to T\ shifted by ±20nK, and have obtained the values 
of dS(t)/d(lnt) about 0.03 and 0.015, respectively. Our final result 0.0233 ± 0.0083 for 
this derivative at t = t* has been obtained by averaging over the three above discussed 
estimates (0.015, 0.025, and 0.03) with the error bars large enough to include all these 
values. According to this, from Eq. (p5|) with £ = 0.6705 we obtain v e ff(t*) = 0.694±0.009 
which, again, is in a good agreement with the value v = 9/13 ~ 0.6923 provided by Eq. (|2l| ) 
at j = 1 and m = 3 and in a disagreement with the RG predictions. 

6 Conclusions 

We have proposed a novel method (Sect. ^|) which allows to predict the exact values 
of critical exponents in the Ginzburg-Landau phase transition model. Our proposal is 
accompanied by a critical analysis of the conventional (perturbative) RG method. In view 
of this analysis (Sect. [2| and ||) and comparison with MC simulation results and experiments 
(Sect. H), our results should not be doubted from the positions of the conventional (RG) 
theory. The best evidence of the correctness of our treatment is the precise agreement 
with the known exact solutions. 
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